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We present a new framework for calculating multi-jet observables through resummation. The 
framework is based on the factorisation of scattering amplitudes in an asymptotic limit. By imposing 
simple constraints on the analytic behaviour of the result when applied away from this limit, we get 
good agreement with the known lowest order perturbative behaviour of the scattering amplitude, and 
predictions for the behaviour to all orders in the perturbative series. As an example of application 
we study predictions for Higgs Boson production through gluon fusion at the LHC in association 
with at least two jets. 



Events with many energetic particles will form the 
backbone of many search strategies at the LHC for 
physics beyond the Standard Model (SM) of particle phy- 
sics, but will simultaneously test our understanding of 
SM processes. Energetic particles of QCD charge will be 
detected as jets, and events with multiple jets require the 
calculation of scattering amplitudes with a high number 
of external quark and gluon (parton) legs. However, the 
calculation of the SM contribution to such processes be- 
yond even the lowest order in perturbation theory is no- 
toriously difficult. Despite recent impressive progress in 
the calculation of many multi-leg scattering amplitudes 
at both tree and one-loop level, the first two radiative 
corrections in the perturbative series (i.e. full next-to- 
next-to-leading accuracy) for observables at the LHC are 
known only for very few cases d, 2, 

One estimate of the effect of higher order correc- 
tions beyond what is currently calculable in full ffxed- 
order perturbation theory can be obtained by interfac- 
ing the fixed perturbative calculation with parton shower 
programs 4, 5], which sum the effect of further soft and 
collinear radiation (i.e. of low invariant mass) to all or- 
ders in perturbation theory. 

As we will demonstrate in this paper, there are impor- 
tant processes where the perturbative corrections from 
hard radiation (i.e. of large invariant mass) is sizeable. 
In such cases, the parton shower approach cannot be ex- 
pected to capture the dominant effect of higher order 
corrections. In this letter we will develop an alternative 
formalism for estimating the all-order perturbative cor- 
rections, which does not rely on the description of only 
soft and collinear radiation. Instead, the formalism will 
utilise results on the factorisation of matrix elements in 
the opposite limit of infinite invariant mass between all 
radiated particles. We will see that without further mod- 
ifications, the results obtained in this limit lead to a poor 
approximation of LHC scattering amplitudes. However, 
known analytic properties of the full perturbative am- 
plitude can be incorporated, such that the approximate 
amplitudes agree well with the full result at low orders 
in perturbation theory. The approximate amplitudes can 



then be applied to any order in the perturbation expan- 
sion, where the full results are incalculable. 

We will validate the approach by considering Higgs bo- 
son production through gluon-gluon fusion (GGF) (medi- 
ated by a top-quark loop) in association with at least two 
jets (the /ijj-channel). We will explicitly compare our 
approximate amplitudes with those obtained using full 
fixed order perturbation theory, where such results are 
known. This process can be used to search for the Higgs 
boson and also potentially to measure its couplings 
to the top quark [7|. The /ij j-channel also has a large 
contribution from weak boson fusion (WBF)[^, [^, [lo| . 
and by suppressing the GGF-contribution it is possible 
to measure the coupling of any Higgs boson candidate to 
the electro- weak bosons, and thus determine 11] whether 
the properties of the candidate match those of the SM 
Higgs boson. The suppression of the GGF-contribution 
is achievedfl^ to some degree by applying specific event 
selection criteria, as discussed in Table HI However, only 
by calculating higher order corrections is it possible to 
estimate the efficiency of such cuts. 

The current state of the art in fixed order perturbation 
theory for Higgs boson production through GGF in asso- 
ciation with at least two jets includes only the first per- 
turbative corrections, and is presented in Ref. |13i]. The 
effects of further soft and collinear radiation were studied 
14| . In the limit of infinite top mass, the coupling 



Ref. 



of the Higgs boson to gluons through a top quark loop 
can be described as a point interaction[l^ |l^, [l3|. This 
approximation was applied in all these studies, and will 
be applied also in the present one, although this is not 
essential to the approach. 

We will start motivating the need for considering the 
hard radiative corrections to higher orders by studying 
the leading order predictions in full QCD for Higgs bo- 
son production in association with hard jets. We apply 
similar event selection cuts to the study of Ref.[i^ as 
detailed in table [T[29] . All of the following results are 
obtained by choosing renormalisation and factorisation 
scales in accordance with the study of Ref. [l4[ , and the 
following values for the Higgs boson mass, vacuum ex- 
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Pj^ > 40 GeV 
\yj.h\ < 4.5 



Vc-Vd < 
\yc-yd\ >4.2 



TABLE I: The cuts used for all results in this letter, in terms 
of rapidities j/i and transverse momenta Pi,±. The suffices c, d 
label cuts that must be satisfied by at least two jets, whereas 
j labels conditions that must be satisfied by all jets; h labels 
the Higgs boson. 



pectation value of the Higgs boson field and top quark 
mass respectively; 

ruH = 120GeV, (0)0 ^ -^,v = 246GeV,TOt = 174GeV. 
v2 

We also include a factor multiplying the effective Higgs 
Boson vertices, accounting for finite top-mass effects |18| : 
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malisation scales as in Ref. [Ij ('301). With this, we find 



(using matrix elements from MADEvent/MADGraph|2lj) 
the tree-level cross-section from the QCD generated hjj- 
channel to be 28ltiiifbi where the uncertainty is ob- 
tained by varying the common factorisation and renor- 
malisation scale by a factor of two. For the three jet 
sample, all jets must satisfy the left-most cuts of Ta- 
ble |T1 but we also require that there exist two jets c, d 
satisfying all the right-most cuts in Table |T1 We then 
find the leading order cross section for the production 
of Higgs Boson plus three jets (hjjj) to be 257l^2o^- 
The requirement of an extra hard jet with an accom- 
panying a^-suppression leads only to a 9% suppression 
compared to the leading order prediction for hjj. The 
Qfs-suppression of the matrix element is compensated by 
the integration over a large phase space for the third jet. 
The large size of the three-jet rate (which obviously de- 
pends on the event selection cuts) was already reported 
in Ref. [l^l , and clearly demonstrates the necessity of con- 
sidering hard multi-parton emissions. 

We will now describe a method for approximating 
the perturbative scattering matrix elements for multi- 
particle production to any order. We start by recalling 
the result of Fadin and collaborators (FKL) 23] , that in 
the limit of infinite invariant mass between all produced 
particles (the Multi-Regge-Kinematic (MRK) limit), the 
leading contribution is given by processes of the form 

n-l 

a{pa) + P{Pb) a{po) + ^ g{p,) + (3{pn) + h{ph) 

i 

where a,l3 € {(7,9,3}, and the partons are ordered ac- 
cording to rapidity in both the initial and final state. 
These processes allow neighbouring particles to be con- 
nected by gluon propagators of momentum (7^, such that 



Pi = li — Qi+i- We have explicitly checked that at leading 
order, partonic configurations which are not captured in 
this framework account for only 0.8fb (0.3%) and 24fb 
(< 10%) of the two and three-jet rate respectively, even 
when there is no requirement of large invariant mass be- 
tween all particles. 

In the MRK limit, the scattering amplitude for the 
remaining configurations factorises, and the results of 
Ref. [23] can straightforwardly be modified to include also 
the production of a Higgs boson. In the case of the Higgs 
boson being produced between the jets (in rapidity) these 
amplitudes take the form: 



n (-^exp[a(qi)(i/i-i -yi)] (iffs C'w(*>'7i+i) ] 



'FKL 



-^exp[&(5i)(j/3 - yh)]CH (%■+!, gh) 

1h 



We choose the A:t-jet algorithm as implemented in Ref. 1 191] 
with i? = 0.6, and the parton distribution functions of 
Ref. [20]. For the strong coupling a^, we choose renor- 



n (-Texp[a(gi)(2/i-i "2/i)] (jffs C'p,('7i,?i+i) 
-!— exp[a(g„+i)(y^j - yt)] (iffs/'"'"+'''"+'ffM,M„+i) : 



where gs is the strong coupling constant (a^ — f'^'"^ 
colour structure constants; j/i, are the rapidities of the 

emitted particles; s — {pa +Pb)^ is the partonic centre of 
mass energy. The factor C^^ is a Lipatov effective vertex 
describing the emission of gluon i. This has the explicit 
form: 



(2) 



where Sai = 2pa -Pi and similarly for Sbi ■ The notation q± 
denotes the projection of a 4-momentum onto its trans- 
verse components. Also, C'h is an effective vertex cou- 
pling the Higgs to off-shell gluons via a top quark loop, 
whose form has been calculated in [23] . The exponential 
factors a(qi) encode the leading virtual corrections (see 
e.g. Ref [2I): 



a(9i) 



9s Nc r(i 



(4^) 
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(3) 



where singularities have been regularised using dimen- 
sional regularisation in = 4 + 2e dimensions, where /i 
is the renormalisation scale and Nc the number of colours. 
The colour factors in equation ([1]) are for incoming glu- 
ons. The form of the amplitude is the same for initial 
state quarks, apart from different colour factors. 

Eq. H]) formally applies in the so-called Multi-Regge- 
kinematic (MRK) limit. Thus, it is clear why this may 
be a good starting point for describing matrix elements 
with many hard partons. Firstly, it can be applied at 
any order in the perturbation expansion. Secondly, it 
does not rely upon soft and coUinear approximations. 
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The multi-gluon emissions have in fact not previously 
been studied directly as implemented in Eq. Instead, 
a simplified version has been used, which is equivalent in 
the MRK limit. In this limit, the squared 4- momenta 
fulfil qf —\q±i\'^, and the squared Lipatov vertices 
-C^^C^ -> 4i^^i^j^l^ii±il^. This means that in the prod- 
ucts of Eq. Il]), only squares of transverse momenta ap- 
pear. Extending these kinematic approximations to all of 
phase space (not just the MRK limit), the sum over j,n 
and the phase space integral over the emitted gluons can 
be approximated by solving the BFKL -eauation\2(y\. In 
this form, the framework has previously been extensively 
applied to other processes. In the present context (after 
implementing local 4-momentum conservation, which is 
strictly sub- leading in the BFKL approach), we find by 
expanding the solution in powers of Ug that the lowest 
order BFKL results for the hjj (554fb) and hjjj (775fb) 
cross sections differ from their full leading order coun- 
terparts by 97% and 200% respectively. The kinematic 
approximations are clearly inadequate in describing am- 
plitudes in general at the LHC. 

Instead, we define a set of amplitudes based on Eq. ^ 
as written, supplemented by the following guidelines: 
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FIG. 1: The 2 and 3 parton cross-sections calculated using the 
known LO matrix elements (solid), and the estimate gained 
from the modified high energy limit (dashed). The uncer- 
tainty band on the LO result corresponds to scale variation 
by a factor of two. 



1. Use of full virtual 4- momenta: Rather than substi- 



tuting qf 



as in the BFKL equation, we 



keep the dependence on the full 4-momcnta of all 
particles. This ensures that outside of the MRK 
limit, the singularity structure of the approximate 
amplitudes coincides with known singularities of 
the full fixed order scattering amplitude. 

2. Positivity of the squared Lipatov vertex: The 
square of the amplitudes in Eq. ([T]) are not pos- 
itive definite, when the effective Lipatov vertex 
is applied to momentum configurations very far 
from the MRK limit. It is here possible to obtain 
— C^^C^' < 0, where the minus sign arises from the 
contraction of the gluon polarisation tensor. We 
choose to remove the contribution from the small 
region of phase space where this happens. 

These modifications combine known analytic behaviour 
of the full scattering amplitudes and the factorised ex- 
pressions obtained in the MRK limit to any order in per- 
turbation theory. 

When these modifications are made, the resulting am- 
plitudes do indeed approximate well the known pertur- 
bative results at low orders, and thus can be reliably used 
at higher orders, where full results are unknown or com- 
putationally unfeasible. Using our approach, we find an 
contribution to hjj of 321fb and a^-contribution to 
hjjj of 217fb, within 16% and 7% of the full results re- 
spectively. In general, we find good agreement in a large 
region of phase space, and the level of accuracy reported 
here does not require any fine tuning of cuts. This is 
summarised in Fig. [1] and one sees that the approximate 
results are well within the uncertainty associated with 
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FIG. 2: Distribution of number of jets, obtained using the 
modified FKL approach, matched to the exact tree level re- 
sults for 2 and 3 partons. The shaded uncertainty bands are 
obtained by scale variations. 



the full results, obtained by varying a common renor- 
malisation and factorisation scale by a factor of two. 

Having validated the approximation at low orders in 
Us (where it can be compared with known fixed order 
results), we now consider results obtained using matrix 
elements with any number of final state partons, which 
at present cannot be calculated using standard pertur- 
bation theory. The divergence in Eq. ([1]) arising when 
any — > is regulated by the divergence of the virtual 
corrections encoded in a. Thus the resulting formalism 
is efficiently implemented in a Monte Carlo generator fol- 
lowing the method for phase space generation outlined in 
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FIG. 3: The transverse momentum spectrum of the Higgs 
boson in the process pp ^ h + njets at tree-level for n — 2, 
n = 3, and the transverse momentum spectrum obtained in 
the resummed approach presented in this paper. 



Ref. [13] ■ Given that one knows the fuU tree level results 
for 2 and 3 partons, however (and they are also computa- 
tionally quick), we have combined these results with the 
approximate matrix elements using a suitable matching 
procedure. We find a total cross section of 499j;307fb. 
The large uncertainty in the total cross section due to 
scale variations however is not reflected in the distribu- 
tion of the number of hard jets as shown in Fig. [2l One 
sees a significant number of events with more than 3 hard 
jets. 

The transverse momentum spectrum of the Higgs bo- 
son when produced in association with at least two hard 
jets is shown in Fig. [3l To the best of our knowledge, 



this is the first report of this quantity, in contrast to the 
completely inclusive Higgs boson pt spectrum, which has 
previously been studied in the literature. The tree level 
2 parton final state predicts a bimodal structure, which 
ultimately arises from the azimuthal correlation between 
the jets. This structure disappears when extra radiation 
is added, giving a qualitatively different behaviour. The 
significant difference between the fixed-order spectra em- 
phasises the importance of considering yet higher order 
corrections. 

In summary, we have outlined a technique, not rely- 
ing on a soft and coUinear approximation, for estimating 
scattering amplitudes with multiple partons to all orders 
in perturbation theory, and demonstrated its application 
to Higgs boson production (via GGF) in association with 
at least two jets. Our technique is based on the FKL 
factorisation formula of Ref. [23], with important mod- 
ifications which ensure that the singularity structure of 
the amplitudes coincides with known all-order analytic 
properties of the perturbation expansion. At low orders 
in as, where the full fixed order result can be obtained, 
our description agrees well, which verifies the trustwor- 
thiness of the approach. It captures both real and virtual 
corrections, and can be applied at any order in in a 
computationally efficient manner. 
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